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We study structure and direct decay of the Giant Monopole Resonance (GMR) using the Time- 
Dependent Energy-Density-Functional method in the linear response regime in a few doubly-magic 
nuclei. In these calculations, a proper treatment of the continuum, through the use of large coordi- 
nate space, allows for a separation between the nucleus and its emitted nucleons. The microscopic 
structure of the GMR is investigated with the decomposition of the strength function into individual 
single-particles quantum numbers. A similar microscopic decomposition of the spectra of emitted 
nucleons by direct decay of the GMR is performed. Shifting every contribution by the initial sinle- 
particle energy allows to reconstruct the GMR strength function. The RPA residual interaction 
couples bound 1-particle 1-hole states to unbound ones, allowing for the total decay of the GMR. 
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I. INTRODUCTION 

Atomic nuclei are known to exhibit a wide range of col- 
lective excitations [l|. Amongst them, giant resonances 
(GRs) are of particular interest in our understanding of 
collective motion in nuclei. Macroscopically, they are as- 
sociated to high energy modes of small amplitude vibra- 
tion of the entire nucleus. Microscopically, they can be 
described, in a first approximation, by coherent superpo- 
sitions of one-particle one-hole (Iplh) excitations 0. 

Giant resonances lie usually at energies above the light 
particle emission threshold. Different kinds of decay are 
commonly distinguished. For instance, direct decay oc- 
curs when one particle is emitted leaving the daughter 
nucleus in a single hole state. However, the decay is 
sequential when two-particle two-hole (2p2h) configura- 
tions are populated by the residual interaction. The di- 
rect decay by nucleon emission is of particular interest 
as the emitted nucleon can give a glimpse on the micro- 
scopic structure of the GR. As such, it is an experimental 
tool of choice to study GRs (see, e.g., Refs. 043 )■ 

Collective motion in nuclei can be investigated with 
microscopic theories based on time-dependent Energy- 
Density-Functional (TDEDF) approaches. In particu- 
lar, the time-dependent Hartree-Fock(-Bogoliubov) equa- 
tions in their linearized version (Q)RPA tS, i7|, have been 
widely used to study various GR strength distributions 
(see [1,0 and references therein). When a proper treat- 
ment of the continuum is accounted for, those methods 
are also capable of describing direct neutron and proton 
decay. 

The goal of the present work is to investigate the 
link between the Iplh structure of giant resonances and 
their direct-decay pattern. The scalar-isoscalar Giant 
Monopole Resonance (GMR) is studied in several closed- 
shell nuclei. The response to a scalar monopole exci- 
tation in the small amplitude regime is computed us- 
ing the time-dependent Hartree-Fock (TDHF) formalism 
[lO| , based on a modern Skyrmc functional |ll| . To study 



GMR direct-decay, the calculations are performed in a 
large spherical grid, allowing for a spatial separation of 
the nucleus from the emitted nucleons to the continuum. 
Similar TDHF calculations have been performed in the 
past with schematic energy density functionals to s tudy 
GMR widths and lifetimes from direct decay [l2|, Il3l |. 
These quantities were extracted from the time evolution 
of the number of emitted nucleons. In the present pa- 
per, new informations are obtained by computing energy 
spectra of the nucleons which are emitted in the direct- 
decay of the GMR. The structure of these spectra is then 
analyzed in terms of single- hole configurations of the A—1 
daughter nucleus. 

Theoretical framework and numerical details are pre- 
sented in section ini Results obtained on several closed- 
shell nuclei are discussed in section Hill before concluding 
in section ITVl 



II. THEORETICAL FRAMEWORK AND 
NUMERICAL DETAILS 

L. Time-dependent Hartree-Fock approach for 
collective vibrations 



The TDHF equation reads 
d 



ifi-^P [h [p] , p] 



(1) 



where p is the one-body density matrix and h [p] is the 
single-particle hamiltonian derived from the energy den- 
sity functional E[p\ using 



hij [p] 



SE[p] 
Spjt 



(2) 



In the TDHF theory, the many-body state is constrained 
to be an independent particle state (e.g., a Slater de- 
terminant of the occupied single-particle states) which 
implies ^ p at all time. 
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The one-body density matrix is used to compute ex- 
pectation value of one-body operators for which TDHF 
is optimized [l^ . Information on vibrational modes such 
as energy can then be obtained through the evolution of 
their corresponding multipole moment. For instance, the 
excitation of a monopole mode | v) on top of the ground 
state induces an oscillation of the expectation value of 
the monopole moment [6, 7\ 



which, in the particular case where Q = F, reduces to 
the strength function. For bound states, the latter reads 



(6) 



while a Lorentzian width appears for finite values of Xi, 
in the continuum. 



F : 



E 



(3) 



C. Monopole vibrations in spherical symmetry 



where al^^ {araq) creates (annihilates) a nucleon in r with 
spin a and isospin q. In this mode, {F){t) oscillates at 
a frequency lOi, ~ [E^, — Eoj/Ti, where Eq is the ground 
state energy and E^, the energy of the excited state. 

The linearization of the TDHF equation leads to the 
RPA 0. As a consequence, TDHF calculations in the 
small amplitude limit contain all the RPA residual inter- 
action which is responsible for the collectivity of giant 
resonances. Using TDHF codes ensures that the same 
EDF is used to determine both the HF ground-state and 
the dynamics thanks to the structure of the TDHF equa- 
tion. It is then equivalent to fully self-consistent RPA 
codes. 

In addition to energies of GRs, TDHF has been used to 
study their stren gth distribution (l5U22j . their life-time 
and escape width jl^ [Tsj , and their anharmonicities (23l - 
[25j . Escape width and Landau damping are accounted 
for in TDHF, however, two-body correlations responsi- 
ble for the spreading width are neglected. Inclusion of 
pairing correlations is possible in the TDHF-Bogoliubov 
theory [1, 0, however its numerical applications [2^ 
are time-consuming and prevent the use of large spatial 
grid as in the present study. 



B. Linear response theory 

The nucleus is assumed to be in its ground state at 
initial time. An external field Vext = ^'n{^)F is applied, 
where e is the intensity of the perturbation and rj{t) is its 
time profile. We choose a delta-function for the latter: 
r}{t) = 5{t). The evolution of an observable 5Q{t) = 
{Q){t) — (Q)o in the linear regime, where {Q)o is the 
expectation value on the ground state |0), can be written 



SQ{t) 



ih 



(4) 



where Q{t) = if t < and 1 if t > 0, {\i^)} are the 
eigenmodes, and c.c. stands for complex conjugated. For 
bound states, Xi, = 0, while unbound states decay with 
a lifetime 1/X^. To analyze the spectral response of 5Q, 
we introduce the quantity 



Rq{uj)^— SQ{t)sm{ujt)dt. 
""e Jo 



(5) 



The monopole moment of Eq. ([3]) is used in the ex- 
ternal field. Assuming spherical symmetry of the ground 
state, only monopole modes are excited. It is then conve- 
nient to decompose the response of the nucleus in terms 
of the usual quantum numbers / and j using the operators 



E 

n,n' ^nij 



{nljm.i\F\jiljmj)alij^^ a^'ij 



(7) 



where the isospin is omitted to simplify the notation. The 
spectral response, Rpi^ (w), defined in Eq. ([5]), can then be 
used to analyze microscopically the strength distribution 
of the GMR. 



D. Spectra of emitted nucleons 

The energy spectra of emitted (free) nucleons can be 
computed by evaluating overlaps between the single- 
particle wave- functions ifnijmj (re) (with degeneracy 2j + 
1 in spherical symmetry) and spherical plane waves 
f^fc/imj (rcr) of energy E = where m is the nucleon 

mass. The latter read 



1 



(8) 



where ji{kr) is a spherical Bessel function for angular 
momentum I and momentum k, ^ij^j (^j 4>) is a spheri- 
cal spinor [23|, and TV is a normalization constant. The 
overlaps between a single-particle wave-function and a 
spherical plane wave with the same l,j and mj thus reads 

cy™, (fc) = E / "^fc'j™. i^'')'Pnljv^, (rcr) dr. (9) 

The probability to find a nucleon with energy E and 
quantum numbers / and j is given by 



^l (^)-E E 



(k) 



(10) 



n mj——j 



E. Numerical details 

The evolution of the occupied single-particle wave- 
function is determined from the TDHF equation in a 
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spherical mesh of radius Rbox = 600 fm, with a radial 
spacing 6r = 0.3 fm and hard boundary conditions. The 
SLy4 parametrization [2S] of the Skyrme EDF [ll| is 
used. In a first step, the HF initial condition is com- 
puted in a i?o = 30 fm box with the same radial spacing. 
At initial time, all single-particle wave-functions vanish 
for r > Rf). 

In principle, the nucleon spectrum should be deter- 
mined at large distances from the nucleus to avoid any 
long-range interaction with the latter. In practice, the 
integral in Eq. ^ is performed from Rq to Rbox- Note 
that the upper limit of the integral in the definition of 
the monopole moment in Eq. ([3]) is chosen to be Rq to 
avoid a divergence of 6F due to emitted nucleons. 

To account for the finite size of the box, only spherical 
plane waves vanishing at Rbox are considered in Eq. 
(Their normalization constant Af is determined assuming 
that they vanish for r > Rbox)- As a consequence, the 
momentum of each partial wave with quantum number 
I can take only discrete values. To avoid spurious peaks 
from this discretization, pjj (E) is convoluted with a nor- 
malized Gaussian distribution in energy with a standard 
deviation of 150 keV. 



III. RESULTS 
A. Detailed analysis of the GMR in ^''O 



proportional to the square of the oscillation amplitude. 
The damping of the oscillation is then a signature for 
GMR decay. Here, the decay occurs by nucleon emission 
as Eqmr is greater than the proton and neutron separa- 
tion thresholds in ^^O. The description of nucleon emis- 
sion is possible if the continuum is properly accounted 
for. This is ensured, here, by the large size of the spheri- 
cal mesh to prevent spurious effects coming from reflected 
flux on the edge of the box [s^ . 

The strength function associated to SF{t) and obtained 
from Eq. © with Q = F is plotted in Fig. [2] (black sohd 
line). The main peak at ~ Eqmb. is surrounded by two 
shoulders at ~ 17 and ~ 31 MeV. In order to get a deeper 
insight into the microscopic origin of these structures, the 
spectral responses of the 5Fij{t) [see Eq. ([7])] following 
the same monopole boost are computed from Eq. ([5]) and 
plotted with colored lines in Fig. [2] We have checked nu- 
merically that the sum of these spectral functions is equal 
to the total strength function Rp, as expected. This de- 
composition shows that the low-energy shoulder and the 
main peak are mainly due to pi/2 and P3/2 orbitals, re- 
spectively. The structure of the high-energy shoulder is 
more complicated. It involves a constructive contribution 
of the Si/2 orbitals while the ones reduce the strength 
in this region. The dominant role of the Si/2 neutron and 
proton contributions in the high-energy shoulder has al- 
ready been noted using similar techniques but other EDF 
parametrizations [ij, 22 1. 



To illustrate the method, let us start with a micro- 
scopic analysis of the GMR direct decay in ^^O. The 
evolution of the expectation value SF{t) = {F){t) — (F)o 
of the observable F given in Eq. ([3]), where the inte- 
gral is performed up to Rq, after a monopole excitation 
in linear regime is shown in Fig. [TJ A damped oscilla- 
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FIG. 1: Evolution of the monopole moment in O after a 
monopole boost in linear regime. 

tion associated to the GMR is observed, with an average 
period Tqj^iji ~ 57.9 fm/c corresponding to an energy 
Eqmr = 2Trh/TcMR — 21.4 MeV. In a coherent state 
picture (see, e.g., [1^), the number of excited phonons is 




10 15 20 25 30 35 40 



E (MeV) 
— Si/2 — Pi/2 — P3/2 — sum 

FIG. 2: (Color online) Strength function of the GMR in ^'^O 
(black solid line) computed from the time evolution of the 
monopole moment 5F{t). The spectral responses 7?f,j asso- 
ciated to the single-particle quantum numbers / and j and 
labelled by their spectroscopic notation are plotted in colored 
solid (dashed) lines for neutrons (protons). 

The next step of our analysis is to study the decay 
properties of the GMR and to investigate their relation- 
ship with the previous microscopic decomposition. The 
spectra of emitted nucleons are computed 2250 fm/c af- 
ter the boost, e.g., when the monopole oscillation is fully 
damped (see Fig. [T]). We use the method described in 
Sec. Ill D1 The proton and neutron spectra for each set of 
I and j are plotted in Fig. [3] We have checked numerically 
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TABLE I: Energ 


ies (in MeV) of the occupied single 


-particle 


states m tne tit 


ground state ot (J. 




s.p. state 


proton 


neutron 




-32.4 


-36.2 


lP3/2 


-17.1 


-20.6 


lPl/2 


-11.2 


-14.5 



that, in the linear regime, these spectra are quadratic in 
the intensity of the boost e [l^l • 

First we note that proton and neutron spectra have 
similar global features, as expected for light N ~ Z nu- 
clei, although the proton spectra are shifted by 3 MeV 
due to the Coulomb repulsion. In addition, the Coulomb 
barrier, which is Be — 1.7 MeV in ^^O, reduces strongly 
the emission of protons at energy E < Be due to the 
exponentially decreasing tunneling probability. We also 
note that there is no relevant emission of Si/2 orbitals. 
In fact, the single-particle energy of the lsi/2 orbitals, 
which are the only occupied Si/2 single-particle states in 
the HF ground-sate of ^^O, are —32.4 MeV for protons 
and —36.2 for neutrons (see table|l]). As a result, the Iplh 
states with a hole in a lsi/2 orbital, and contributing to 
the GMR microscopic structure, are bound for neutrons 
and below Be for protons. Indeed, the different struc- 
tures in the GMR spectrum (see Fig. (2) are located at 
energies below ^ 33 MeV, which is not sufficient to emit 
a lsi/2 single-particle. 




E (MeV) 



— Si/2 — Pi/2 — P3/2 — sum 



FIG. 3: (Color online) Neutron (solid lines) and proton 
(dashed lines) direct-decay spectra of the GMR in ^^O for 
each set of single-particle quantum numbers I and j, given in 
spectroscopic notation. 

Let us now attempt to reconstruct the GMR strength 
function from the spectra of emitted nucleons shown in 
Fig. [31 In a pure harmonic picture, the GMR is a co- 
herent sum of Iplh states with a difference ~ Egmr 
between the particle and the hole energies. However, the 
nucleon emitted by direct decay have only the energy of 
the particle-state, which is obviously smaller than Eqmr 
because the energy of the hole-state is negative. To re- 



construct the GMR strength function, it is then neces- 
sary to shift the energy of the emitted nucleons by the 
energy of their associated hole-state, that is, the energy 
of the initially occupied single-particle state. The latter 
are given in tableUfor ^^Q at the HF level. The "shifted" 
spectra of emitted nucleons for each set of / and j, noted 
plj'^\ obtained from this procedure, are plotted in Fig. 01- 
b, together with their sum. It is striking to see that the 
shape of the latter is very similar to the strength func- 
tion shown in Fig. |21 and recalled in Fig. Hl-a, although 
the two spectra have been obtained from totally different 
quantities, i.e., the spectra of emitted nucleons and the 
time evolution of the monopole moment. 

Similarly to the microscopic decomposition of the 
strength function obtained from the time evolution of 
the monopole moment (see Fig. [5]), we see in Fig. Hl-b 
that the low energy shoulder and the main peak are as- 
sociated to pi/2 and P3/2 orbitals, respectively. Again, 
the situation for the high energy shoulder is more com- 
plicated. Although it is mainly due to si/2 orbitals in the 
monopole response (see Fig. [5]), it is in fact associated to 
the emission of p-orbitals in Fig.|4l-b. As discussed before, 
this is because the energy of the GMR is not sufficient to 
promote a lsi/2 single-particle to the continuum. How- 
ever, the RPA residual interaction couples the responses 
associated to the different single-particle quantum num- 
bers to produce a collective monopole vibration. As a 
result, the vibration of s-orbitals can be coupled to p- 
unbound-states contributing to the direct decay. Such 
coupling might be the reason for the complex competi- 
tion between constructive and destructive contributions 
in the high-energy shoulder observed in Fig. [21 



B. Tin isotopic chain 

Similar analyses have been performed for the GMR in 
ioo,i32gj^ nuclei. The strength functions obtained from 
the time evolution of the monopole moment are plotted 
with black solid lines in the upper panels of Fig. |S1 and 
|6l The reconstructed GMR strength, obtained from the 
emission spectra are plotted with black solid lines in the 
lower panels. As in the ^^O case, the agreement between 
the two methods is excellent. 

The gross feature of the strength is similar for the 
two isotopes, with a single peak centered at Egmr = 
17.2 MeV for loogn and 15.8 MeV for "Zgn (see ta- 
ble |ll|, although their widths vary more importantly. It 
is striking, however, that the decompositions in terms of 
single-particle quantum numbers exhibit some important 
differences between the two techniques. Although both 
proton and neutron single-particle orbitals participate to 
the vibration (see upper panels in Fig. |5land|6]), the de- 
cay only occurs by neutron emission for ^'^^Sn, whereas 
^"^Sn decays through proton emission only. This can be 
understood by the fact that the proton (neutron) sepa- 
ration energy increases (decreases) with the number of 
neutrons. These quantities are given in table |TT1 together 
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FIG. 4: (Color online) (a) Strength function of the GMR in 
^^O obtained from the time evolution of the monopole mo- 
ment, (b) Spectra of emitted neutrons (colored solid lines) 
and protons (colored dashed lines) "shifted" by the energy of 
the initially occupied single-particle state (see text). Their 
sum is shown in black solid line. 



FIG. 5: (Color online) (a) Strength function Rf (black solid 
lines) of ^""Sn. Its single-particle decomposition, multiplied 
by a factor of three (-Rf,j x 3), is represented with colored 
solid (dashed) lines for neutrons (protons), (b) Spectra of 
emitted particles "shifted" from their single-particle energies 
are plotted in colored solid (dashed) lines for neutrons (pro- 
tons). Their sum appears in black. 



TABLE II: Neutron and proton separation energies and 
Coulomb barriers in Sn isotopes. The GMR energy corre- 
sponds to the position of the peak in the strength function. 
All energies are in MeV. 





Sp 


Be 


Sp + Be 


Sn 


Egmr 


^"•^Sn 


3.1 


8.1 


11.4 


16.9 


17.2 


"2Sn 


15.6 


7.4 


23.0 


7.7 


15.8 



with the Coulomb barrier for protons and the position of 
the peak energy associated to the GMR energy. Indeed, 
the binding energy of neutrons in ^*'°Sn is of the order 
of Egmr and proton emission is favored, while the com- 
bined effects of an increasing Sp and the Coulomb barrier 
Be prevent proton emission in ^'^^Sn. 

As in the case of the high-energy shoulder in the ^^O 
spectra of figures [2] and |4l the difference in the micro- 
scopic decomposition of the GMR strength with the two 
techniques is a signature of the RPA residual interac- 
tion which couples the bound particle-hole states to the 
unbound ones, allowing their decay. In particular, the vi- 
bration of the protons (neutrons) in i32Sn (looSn) decays 
via neutron (proton) emission thanks to the collective 
vibration of both proton and neutron mean-fields. 



T ' ' ' I ' ' ' r 
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E (McV) 

— sum Si/2 — Pi/2 — P3/2 d3/2 dg/a 

h/2 h/2 §7/2 §9/2 hii/2 

FIG. 6: (Color online) Same as Fig.[5]for "^Sn. 
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IV. CONCLUSIONS 

Using a simple single-particle decomposition, the mi- 
croscopic structure and (direct) decay of the GMR in 
some spherical nuclei has been studied. The GMR 
strength was obtained by two techniques: (i) from the 
real-time small amplitude monopole response and (m) 
from the spectra of emitted nucleons. In (ii), each micro- 
scopic contribution has to be shifted by the single-particle 
energy of the hole state. 

Although both techniques lead to almost identical 
monopole strengths, their microscopic structure (i.e., 
their decomposition onto the single-particle quantum 
numbers / and j) may be very different. This is under- 
stood by the fact that, although GMR energies are usu- 
ally above the nucleon emission thresholds, only part of 
the particle-hole states lie in the continuum. Those which 
are bound are coupled to unbound states via the RPA 
residual interaction which is responsible for the collectiv- 
ity of the vibration. This is illustrated, e.g., in monopole 
spectra of neutron-rich tin isotopes: although both pro- 
tons and neutrons contribute to the monopole vibration, 
only neutrons can be emitted because protons are more 



bound and have to overcome the Coulomb barrier. 

A major improvement of the method would be to in- 
clude two-body correlations, like pairing interaction and 
collision terms. The latter are known to be responsible 
for the additional spreading width of giant resonances as 



well as a fragmentation of their strength function [31 



They are also expected to modify the structure of the 
spectra of emitted nucleons by coupling Iplh states to 
2p2h. A first step would be to use the Extended- TDHF 
[Hi-Isl or the time-dependent Density-Matrix 
formalisms. Because of computational limitations, spher- 
ical symmetry (and thus, only monopole vibrations) 
might be first considered. 
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